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ABETHACT 


The object of this thesis is threefold, to present 
the theoretical anslogy betwee cn the floes of a comeressible 
g&s, such as air, enc the flew of water with a free surface, 
to deseribe the détails of construction of the seter table 
hich utilizes the above enalogy, and to outline ané demon- 
strate the applications of the water table in the investi- 
gation of air flow, 

The anmilogy between air flow end water sith a free 
surface is presentec in Part I. Part II contseins a ceserip- 
tion in detail of the construction of the water table at 
the Rensselsxer Polytechnic Institute. The appiications of 
tration of its use are found 





the weter table ond one demon. 
in Part III. 

This work wes done during the months of Februarg 
through Mey 12947 at the Rensselaer Polytechnic Institute. 








ROR INCLs TURE 


The fcllowing nomeneleture is ised in this pepper: 
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5 speeific heat «t constant pressuré, Btu pwr ib per deg F 

Cy + specific heat at constent volume, Btu per lb per deg F 

‘se @ = totel differentiel 

cme tee 

= 773.88 ft-lb per Btu 

te sechenioel work, Bou per 1b 
ee 

| Ek Woe poe be 

1 = gut constents 98.8 for aie = 

static temperature, deg Rankine 

eee 

= velocity, foe 

: = velocity of # sound wove, fps 

¢ maximum velocity 

‘Wz Weight flon, 1b per see 

Ua Pato oF gus apeeksie nests, 6/0, 5 1.969 for vole ott 

Po mse censity, slugs yer ou ft 





















































The following nomenclature is weed in this papers 


ee wont of the veloeity in the x,y, 4ifeetions 
a Water depth, ft 

. A total head (water cspth when v=0), ft 

=: , gatel Kea after hydraulic jum, ft 


._ 2 

















pressures, pafa 
te = | quantity ef flox, eu ft per see 
velocity, fps . — 
“welecity of prepegstion of & surfeee exve, tps 
reetangular coordinates in flow space 
yP yet, 2 , state (subseript *o*) 
(- mast @ensity, slugs per ev ft ’ 
es 
Pn. ~_ 2a 
=—— : 
ae 
°c a . = = 














LRTHRODUCTICR 

The operation of « supersonie Wind tummel for test 
work at transsonic and supersonle air speeds is an expensive 
prajeet, re ulring costiy e:ulpment, elaborate instruménte-~ 
tion, and en enormous amount of power. A elzeable steff of 
skilled technicians and much time are neeéed in the fabrica- 
tion of models to be used in the tummel and in the «e¢tual 
operation and maintenence of the turmel ena its edjunets. 

The reduction in the emount of test work thet must be 
done in high velocity air would save « considerable emount 
of time and expense. This would facilitate Cesign, @xperi- 
mentation, and theoretical research in the fleid of high 
Velocity airflow, 

It is the intent of this paper to show that the water 
table miy be sdventsgeouSly employed im conjumetion sith 
the Wind tummel. Its use for preliminery work would replece, 





to a2 considérable extent, the oyeratioans 
wind tunnel. 

The wark herein vresented wes Wwerteken With three 
primary objectives in mirc: (1) em investigation of the 
eListing theoretical] bnalogy between the two dimensional 
flow of a compressible gas, such as air, end the flon of 
shallow water with e free surface, (Z) the construction 
of & Water table, & devic® Designed to utilize the flow 
enelogy, ami (3) an investigation of the flow sround « 


necestery in the 


fee basic models in orcer to check the anelogy *tith 
existing eir flow theory ard to demonstrate the usefulness 
of the Water table, 











By # compurison of the energy emt continuity equations 
for two dimemsionel gas flow and water flow with « free sur= 
face, it was found that an «nalogy exists between the ratio 
of absolute air temperetures end water desth ratio, the 
ratio of gas cemsities and water depth ratio, and the ratio 
of gas pressures end the syvarse of the ®ater depth ratie. 
This analogy ie aquentitatevly correct for « fictitious gas 
having e# ~ C67, of 2.0. "Hydreulie jump* in water flow 
is used to represent “eoosression shock" im sir flow, 

Im order to cemonstrate the usefulness of the weter 
table as epplied te the study of supersonic eir flow, it 
was originally intended to check a few of the existing 


_ #rimeiples of supersomie air flow theory on the weter table. 


Experiments were siso proposed to investigate that part of 


high velocity air flea ebout which Little is known «at 


aterecd in the construction 





present. Due to diffieulties encou 
of the water table, littie time was available for model 
sonstruction and test work. Therefore, oniy one spplication 
of the water table is presented. This is am investigation 
of high velocity flow urow a fifteen degree wedge. 
This work was cone cguring the months of Februsry 
ough Bey 1947 et the Rensselaer Polytechnic Institute. 





fhe water table was constructed in the Bechenical Sngineer- 
ing Leboretory cf the Pensselaer Polytechnic Institute. 








rae L 
THe #aATER ANALOGY OF SUPRAOORIC ALP FLOW 


Reversibility, in the case of & someréssible fluid 
suoh as air, is defined by the stetement thet the chonge 
im internal energy (Cy@T) is caused entirely by exsension 
or compression work (APa¥) or 

S(CyEP) + PAV = 0 (1) 
fhe state of s eompressidlie fluid is defined by the gas 
equs tion: 

PY = RT (2) 
If the gas equation (4) is differentiated we hever 

POV+VaP = RET 
and since, at amy instent, V= Of, then 


Pa + Rage = Rat 


pey 4 ROT - Aree (5) 
Substituting this value fer PdV in evuation (1): 


(C+ Ser 7 arse (4) 
but 
= C, += C, = AB mere me 


fhen esuation (4) becomes; 


. oe ot = 1 


fan © 
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Integrating (5) gives: 


or 

7 v7 
o~ = Go) (8) 
From the gas eyuation, using (#) pownds of aire 


PY = PRT 


ary, 
~3 
Seat 


? = P ext 
Using this P in eyuation (6) gives: 


it 


L-@o (8) 
By combining a and (8)x 
eG @ 
Equations (6), (8), end (G) define the relationships 
betweer P,f/, and fT for compressible gas flow whieh is 
reversible (no hest added — no friction). 
The general energy ecuation for gas flow: 
AQ + dL = Oye +.GPd¥ + ANG? +hvéy 
For the case of no heat added and no outside work done on 
the ges (dQ ~ 0, dL = 0): 


CyGT + WGheV + SVOP + sovdy = 0 


Then: 


but 
| t Pay WoP « oT 
Cyt + ipa + NAP = Cy 
~ 
ja - ~ Lé 
Gpat + sovdv O (10) 








Integrating between the Limits of T, and T: 
a id PS 
If fp is taken et the state of rest (¥, = 0), equation 


(10) becomes: 


Lg 
Cp(T> -T) = ZygY 


Ueing Cy in the units of foot-sounds/.egree Rankine: 
C)(T>) -T) = Bev or v's 2g0,(t) -7) az) 
The energy ¢.vation in the case of frictionless 
eater flow states that the sum of the potential energy 
and the kinetic energy of s weter partiele is constent, 
Consider a flow filament (Fig. 1) which passes 
through the point y,, 2, of the initial cross section 
x - Oe Along this filement, between the pressure P and 
the velocity v, the energy a vvation iss 
P + lev +egz = Constent - Pot bev .tee2, (28) 





On the free surfsee of the zater P is constant and ecuel 
to the atmospheric sressure, Pa. In whet follows, we may, 
without error, set this to sero since only differences 

in pressure cre of physies] significance in the case of 
incompressible flows, If the weter flows frem an infinitely 
Wide besin then v, = 0. The curvature of the free surface 


ose this 





is also zero at this point. It is logieel to Fre 
point 2s the reference point x,, ¥,5 %,- The corresyonding 


water depth is denoted by h,. For the sbove reference 


point, the Bernoulli e..untion is: 
P+ hev'+ e&2 = Po + PB%, 









7 ; 23 » s 





«<a cs = S&S 





or 
v = Re{z, -2) - a, Pi/ (14) 
It if reasoneble to Sxswme that the verticel 
aeceleration of the water is negligible compered with 
the gccelerstion of grevity. Wider this assumption 
the static pressure at « point in the field ef flow weries 
linearly with the vertices] distence from thet point te 


the free surfx@e: 


P= pe (h,~z,,) (35) 
and 

P = pg(h-s) (16) 
Substituting (15) end (16) im (14) gives: 

¥ = 2g(h, -b) = 2gah (27) 


The above energy equetion is valid fo: the flan 
filament vassing through ¥, end 2, at K = 0. Binee, at 
x =O, all the filaments which lie one ehove the other 
have the seme h, and v¥, (zero), end since equation (17) 
does mot contain z, the velocity v at x and y is constant 





ever the entire depth and is e.uel to the difference in 
height between the total head (ho) snd the free level 
(h), (ah), at most cen eyuel h,. The maximum velocity, 


thercfore is 
Vnax = | tah, 


If we write ejuction (17) in dimensionless form: 
aWy* ~ ABR 3. 
Fae So - i = 














From @yjuation (12) it is obvious that the meximum 


velocity im © gas iss 


end (19) 





* | : 9 

From eqtatibes rey silt ( (18) 3 it may be Seen thet 

the ratio of the velocity to the maximwm velocity for 
the water and the gas flo#s becomes evel if 


L-grel-<¢ 


“a2 ") 
Bquetion (0) shows thet, with respect to the 
velocity, there exists en anelogy betkeen the two flows 
if the depth ratios are compar®d fith ebsolute gas 
temperature rotios, 
| By comparing the eyuation of continuity for the 
®eter flow with a free surface to thet of a gas, it is 
possible to obtcin smother enelogy. 
. rat x,y, © smell fluid prism with the 
Gimensions du, Gy, «md h (Fig. 2). Let a, b, and ¢ 
represent the components of the velocity in the direction 
ef the x, y, end 2 Axes respectively. 

If the assumption is made that the vertice1 eccelera- 
tion of the water is negligible in comparison to the 
aeceleration of grevity, evuation (16) cen be used in dif- 
ferential form os: 














2 .,,.d8 
ox “PE x 


and 

by =C8ay 
The tight sices of the above equetions sre independent of 
2, Which mecrw that the horizontal asaeeelerstions of all 
points alomg a vertital aré sliso imciependent of sg. The 
horizontel components, e end b, are then constant over the 
depth h. 

The continulty eaustion for this type of weter flow 
states thet the rate of mass flow into the prism is exsual 
to the rate of flow out of the grism., Ginee the danmsity of 
the water is censtunt, the infloring volume per unit time 
(49, _.) must equal the outfloring volume (do 

co, = ehdy + bhdx 


out? 
In 
ead | 

CQaug = (@ MMEax) (besBsxday + (oa ghay) (ne Shay )ax 
By #xXpancing and neglecting infinitely small magnitudes 
of higher order, we have 





Ou CEAy + 2\tb) = 0 
- | 
ant dividing by dxdy: 
A (Mb) «6 2 
Dihe) , Mab) « 6 (22) 


Equation (21) is the sontinuity equation for 
Stationary weter fics, 

Tre eontinuity e.uation for a two cdimensicnel 
compressible gas flow iss 











{-i5 


oo + ee « © (xe) 


Equations (£1) and (se) obviously hewe the same 


a? 


“ta 


form. From these @)uations it is possible to @erive « 
further condition for the snalogy of the two Flows, 
namely, thet the aensity of the ges flo corresponds 
tc the depth, h, of the water, Expressing this ema logy 


aS a Gimersionlgss ratie;: 
Ro Be 
it is now .,cesible to investigate the physiésl 


7 

ne 

ui 
ae” 


mature of this gas which we «re comparing with the flow 


of water with « free surface, The sciabsati¢ aeyuetion (8) 


Ea? 


From (20) and (23) we have: 
—_—— _ -. a 
By substituting these values im equation (6), we have 


2)" 


whieh obviously is satisfied only if 
Thus it is shorn that the flox of the weter is quamti- 
batively comparable with the flow of « gas having « retic 
Ja She qual to #0 
Amother ansiogy, probebly the most important from 
an experimental pcint of vies, is obtained from ine ga8 





states thet: 


the e.jmwation: 





equation: 








or 


i re 
te = €% 
Le > 


Ny substituting tha val of 2 ame a im beri of 
@ 


Co ° 


fe Raves 
= gh (wa) 
Roumtion (84) cen alec be ¢Cerived by using olvwer 
of the eGiabstic eyuticens (6) or (3) ond Y= 20 
Por 2 compressible gas, the velocity of propagation 
of 4 sound weve v2 * - |/raht.1 the velocity of groyage 
tion of a eatune way tn wheiipe sep fe i = =V me 
mss velocity 1s culled theveriticsl velveity*. 
In © gas expending through ® nostle there is a 
ritiesl] pressure retic «t mhich senie velecity occurs = 
(St the throst, In the case of sir with» J= 1,385, — 
Se ree ee one eae 
See Prsees opie euatkn Gbrngs, Se Pi 
0 ) is less then this value, sonic velocity 
yt occur in the nozale threst. Am enelegows condi- 
. Occurs in the water Plows. ier 
| Mme erstieel velocty of sever ds eb. Prom the 
By equation (17) wo have: 
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a Ai 





“Held p, betiey, "The therwosymanton of et sigh 


*gonliention of the ¥ethods of Gas 
“BES metoctt ath cia peee pores mt 
O84, p. 8, 1840, ; 

















vo « eg(h, ~h) 
By substituting v,, = V gh, eiuvetion (17) becomes 
gh - Ugh, - “gh 
or 
4a 
Them at amy voint where the weter denth is tee thirds of the 
total heac, the water is flowing #t the orliiesl veloeity. 

In air, if the #ach musber e: >) is les? than one, 
the flow is seid te be Msubsonie" and if His greater then 
one, the flow is *supersonie"™, In water, ir is less 
‘them one, the flox is kneen as "streaming flow end if + 
is greater then one, the flow is "shooting flo#*, 

If air is flowing ot supersoni@ velocity, a disconti- 
nudity esy oceur in which the velocity suddenly decreeses to 
6 subsonic value which setisfics the conditions of flcw,® 
This discontinuity iz eslled » "plane compression shock" 
and theorwtieally takes place in « very short distance, Just 
as the velocity suideniy decreases through «a plane shock, 
the pressure sudcenly incresses, then supersonic gas flor 
is foreed tc chenge its flow direction due to an obstruction, 
an ®angle shock" will oceur.* it is Known that am enalogeus 
discontinuity, called the “hydraulle jump*®, mey oeeur in 
Shooting weter flow. Asin a gas, two cases of the jump are 
possible: (a) In the right Jump, shcoting water is con- 





3. Weil P. Batley, "Abrust Mmergy Trensformations in 
Flowing Geses*®, appendix B, p. &. 
4. ltbié, Ds ots 
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verted into streaming flow and (b) im the slant or oblique 


jume, the flow mey or a mot go to etreaming efter the 





jump, depending on the a 


i 


of the water, 


& swemarg of the flow analogy is as follows: 


THO DIRSBCLORLL 





8 FLOW 





flew mediun With X= Eyd 
boundsries similar 


















Velocity 5" 


Temperature FEtLO gg 

o 

Density retio,f 
Qo 

Pressure ratio, gf 

Q 

Sound velocity, ve 
Hach Number, > 

fubsonic flow 


Supersonic flow 





Compression shock 
(vlane and angle) 


LI.UID FLOW RITE FREE SURF, 
I! GRAVITY FIeLp 


incompressible fluid (wate. 


fide boundary vertical 
Bottom horizontal 


Velocity = Poee 
mix OF 


Depth ratio, 


ae 


Depth réetio, 


S@uere of acy va Petios (py 
@ave Velocity, v¥ 

Nach number, + 
Streaming flow 
Shooting flow 


Hydraukic Jump 
(normal and elent) 











Thus far we have shown the amalogy that can be 
made between water end s ges having s JVof a0. The analegy 
between the compression shock and the hydraulie jump, 
however, does not strietly held. The energy eqguaticn (17) 
between the veloelty and the depth for suter flow ise 

. «= eg (My —h) 

where the total head, (ho }, is constant. Im the case of 
the hydraulic jump a portion of the xinetic energy of the 
water is Bomverted into heat. For this reason the tetal 
head efter the jump, (hd), is smaller then the total head 
before the jump, (h, ). After the jump the energy equetion 
ise 

- “ = hg (ha-h) 

The energy loss during the jump beers s timple relation 
to the intensity of the jausp. In the flow over « horizen- 
tal bottom the potential energy is a minimum if the water 
Gepth, kh, is vero. For 2 mezs of witer, m, at 2 depth, hy, 
the potentizl energy is P = mgh/2. 

@inee the kinetic energy st points of rest is e2tal 
to £ere, the energy lose, (4e), ehich oeewrs in the hycreulie 
jump may be computed as the difference of the potential 
energy «t a point of zero velocity before and after the 
jump, or 

Ae - mg (Bo ~ a (25) 
By dividing this energy loss by the energy before the 
jump, @ = mane, the reletive energy logs is obteined #5: 
Es BE w I~ hi/he (26) 









<< = _ 
~~ — er: “ 
a) ese a. ole ee Oe 
= — 2 5 
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This is the relative amount of energy which hes been 
eanverted into heat and is lost energy insofar 88 water 
is concerned, 

In 4 gas, the heat genereted by a shock is not lost, 
but is merely converted into thermal energy, end the total 
temperature, consequently the total energy, i8 the seme 
before ami sfter the shock, Singe the eomperable weter 
magnitude to the gas tempereture is the denth, the analogy 
is not strictly true after a hydraulic jump. The energy 
Loss in water? is extrem@ly small over « lerge region of 
jump intensity, the relative loss being less than ohne 
pare tee ee te 3.0. As &© rewult of this small sheek 
energy less, the anslogy of the two types of flow ie still 
Walid as a very ¢lose sporoximation within the range of 
Bech numbirs eurrently employed, 

& Preiswerk, Op, cit, Part II, p. #1-He 
t= 





PART II 
CORSTRUCTION OF THE E.P.I. SATHR TABLE 


The essential element of a water table is a horizontal 
surface with vertical sides over which « shallow stream of 
water flows. The necessary edjumcts to the surface are e 
framework to support it, the tanks, pump, and piping to 
supply the flowing water, and the necessary measuring instrv- 
ments, The follcewing description of the water table construc- 
ted in the Mechanical Fngineering Laboratory «t the 
Rensselaer Polytechnic Institute will serve to illustrate 
one method of fabrication, 

The surface is a sheet of plate glass, $" x 48" x 60", 
Glass was selected for two reasons. First, glass is one of 
the most frictionless surfaces obtsinable sufficiently rigid 
‘to maintain a horizontel plane. Second, the use of glass 
makes possible the taking of photographs of the flow either 
from above or underneath the surface of the table, illumina- 
tion being supplied on the side opposite to the camera, 

The frame is constructed of white pine. The table box 
is constructed of £" x 6" plenk, the legs of 4" x 4" timbers, 
ané the crossbars end leg braces of «" x 4" stock, An 
isometric drawing of the frame is shown in Fig. 4, and 
actual dimensions cen be found in Fig. 4. All joints ere 
mortised and securdd by bolts, To provide for leveling the 
table, the bottom of each leg is fitted with a 2" steel 
plate imto fwhich is threeded a @" steel bolt. 


3089 
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The giass is suppertec on the sides by L* x 2m strips 
of wooGd Screwed to the frame omi the forserd amd sfter enis 
rest im grooves cut in the «" x 4" erossbers, Beteils of the 
support arrangement are shown in Fig. 7. Two 1" angle irons, 


surfaced with masonite, placed one-third the Length of the 





giass from each 
glass is held firmly in pleee by 1” x 2" wood strips run 
timc length of the table, These strips ere sovered with 


end, supply Support across the table, The 





cepper Tinshing in erder to prowide smooth vertieal walis 
ana to prevent the sozsibility of warping due to mwster-to= 
woes contact, To imsure wetertightness emi cvoid » metal- 
te-glass contact, o strip of i" sponge rubber is secured 
between the copper and the glass, 
- fhe tanks are conntructed of lisgauge rolled sheet 
steel. Dimensions of the forward tank are shown in Fig. 6. 
This tank has a capacity, up to the lip, of about five cubic 
feet, The forwerd tank is bolted to the freme om three sides, 
fhe four-inch lip rests upon the forward &" 4 €* erossber 
emda prejeets over the forverc eage of the glass. The iip is 
prevented from direct contact with the glass by a Strip of 
4* sponge rubber, 4 three-ineh flenge with « standard 1g* 
pipe thread is braseé to the bottom of the tamk to «ecom- 
modete the sugply pipeline. 

Dimerisioms of the ofter tank are shown in Fig. 6. 
It has « capaedty, up to the lip, of about 3.3 evbie feet, 
The after tank is bolted to tne after end of the frame and 
the ome-ineh lip rests im a groove cut im ithe efter :* x &* 






Meter eee - 


aie "ee, -= 
time > ~~ width armre )* 


. 





> a 






if : 
7 a 
e = > ——-?_> 





crossbar, A strip of co; per flashing runs from the after 
edge of the glass, over the tamklip, ami into the tank. This 
provides a smooth surface for the water flow amd prevents 
water-to-sood centact, Three four-inch flanges with stendard 
two-ineh pipe thread are brased to the bottom of the tank to 
accommodate the dischirge pipelines, 

The water for this table is self-conteined and recir- 
ewlating. For this pursoge there is an SO-gellon capacity 
sump tenk ct the after end of the teble. The pining and 
eomnections between the sump, two tuble tamka, and sump tank 
are shown in Fig. @. Sinee tne pump is an elternmeting cur- 
rent electric-criven constant speed centrifugal tyre, a 
eomtrolling vwilve on the discherge side of the pump and e 
Mater by~pass line are used to throttle and divert the 
flow, All piping is two~ineh steel exeent for the by-pass 
line ana short section lesding into the forvard tank, which 

are 1¢" steel. Valves are loested im the three coischarge 
Limes from the efter tank to regulete beck presturé on the 
table surfece, 

To maieure water flow, there is « flat plate orifice, 
1.25" in diameter, in the ig" line leeding to the formaré 
temk. Pressure tacs, locatec in aeccordanee with the £SRE 
research publication on fluid meters®, are connected by 
rubber tubing to a U-tube weter-filled menometor. The 
orifiee and manometer Were calibre>ted by weighing the 
are SRE TER Tee emntinstieener 


& "Flinid Meters, Their Theory and ApOhlieation,*® Part I, 
ASWE Reswareh Publication, LfS1, 
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amount of water flowing in e measured intervel of time, 
A cellbration curve is shown im Fig. Gs 

To measure depth om the water table surface, a hock 
gauge, capable of measuring to 1/1000 foot is used. The 
hook Was replaced by e straight pointed rod for grester 
convenience and sccuracy. The gauge is mounted on e wooden 
block which slides freely alerg e 8" x 2” wooden stringer, 
The stringer rests om the teble frame, Depth is measured by 
teking the cifference of s reading at the glass and one at 
the weter surfsce, 

A wire basket filled eith smali rocks is sed to 
smooth out the water flow, This basket is fabriested of 
two stripe of 4" wire mesh which are three inches apart. 
The basket rests im the center of the forwerd tamk and 
runs the width of the tank. It also serves the purpose of 
& water filter, Two weirs sre provided te supply the neess~ 
sary totel head. They are fabriected of 14 gsuge shevt 
bress, rolled an! bent to the desired shape. The uncernesth 
Side of each weir is fitted with three equelliy spaced 
‘Wooden bloeks into each of which is seeured a 3" bolt. The 
Wooden bloeks arw held in glece by smell bress sechine bolts 
“whieh aré countersunk into the ten of the weir, The bolt 
heads are covered with solder which was smoothed to conform 
to the shape of the surface of the weir, The securing bolts 
projeet through holes in the lip of the forward tank and 
erostbar anc. are secured by nuts besring eguinst the under-~ 
sige of the crossbar. Natertightness is mainteineé by using 
eaulking compound between the weir and the teble sides and 
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tamk, Details are shown in Fig. 7. 
| Tre wooden acdels usec thus far in comgumetien sith 
experiments on the table were fabriested of white pine end 
finished with sper varnish, semded to a fine polish, 

To minimize surisce tension, a wetting egent is used 
in the Water. 

Fig. 10 is « photegreoh of the completed table in 
whieh the frame, pump, piping, valves, manometer, uné éapth 
gauge are visibie, 

Fig. 11 is e photograph of the table tep leoking for- 
werd, in whieh the gless, weir, anc rock filter are visible, 
A typleeal experimental set-up is shown in plece om the 
glass. 

- In the course of construction bwe major diffilcuties 
sere encountered, mounting the glast end obtaining laminar 
flor. 


. The originel intention ms te use #" plate glase both 





fer rigidity and strength, This gless soulé not be obtained 
within the time evallabie. After investiestion of strength 
and rigidity properties, §" olate glsssa vas determined te 
Be Satisfactory ami nes instelled. This glass was supported 
only on the four sities. tbout ome month after instellutien, 
this glass fractured, Apomrently, the glass had sagged 
Slightiy in the center and the stress set wp, in time, re- 
suited in fracture, A 





ond 3" plate glass was instelled 
and the two angie iron braees were added ws edtitional 
support. after ebout three seeks, 2 smell fracture S468 


experienced at the after corner of the table. The c#use of 
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this seeomi breek ig not definitely known. No impact spot 
is discernible anc the glass is perfectly bevel, This break 
does not affect the operation of the table. a sheet of 3° 
piate glass hes been ordered and will be inetelled zs soon 
ag delivered, 

The Reynolds number of the flew on the tabae is very 
near the criticel. In order to obtain Luminar flow, the 
water must proceed from its entrenee inte the forward tank, 
through the rock filter, over the weir ané om to the table 
with gredueal turns, The selr surfece must be extremely 
Smooth since any proetuberan@e #111 set up a disturbence 
end foree the water into turbulent flow, The firet weir used 
as comatructed cf weoé. A reslly smooth surface coulda rot 
be obteined on the wood available, anc this wes disearced 
in favor of the weir constructed of brass. Leminar flow is 
ties. 
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out the vracticel renge of flow veloei-+- 














PART Il 
APPLICATIONS OF THE WATHR TABLE 


The anpplicetions of the water table in connection 
with the study cf air flor ere many in number. The water 
table may be used to study any type of tro-—dimensicnal 
air flow which Goes not involve a change in the total 
energy. Some of the more important applications, in the 
opinoin of the *riters, are as follows: 

(1) the investigation of the flow of supersonic air 
in and eroumé cirfoils, inlets, nozzies, and diffusers, 
Modeis must be designe@ smi fabrieceted for air flow with 
a Ye 2.0. 

(2) eonfirmation of the validity of assumptions mead 
in the development of air flow theory. Using Ya %.0 in 
the e uations, the theory can be checked with experimental 
results obtained from the *ater tabie. The water table 
offers en excellent medium for studying the little known 
transsonie region, 

(3) as a laboratory apparatus for demonstrating the 
accepted theory of various types of gas flow. 

One exemple of the appliestions of the eater table 
is demonstrated in the case of flow of supersonic eir 
#rouwnd a simple eecdge, The surpose of the experiment wes 
to check the ester enalogy with the existing theory on 
angle shocks .? 





7. Bailey, Op. cit, p. 4. 














A sketeh of the wedge used in the exberiment is shown in 
Fig. & Experimental data can be found ip Table I. Table ITI 
eontains the results caleuluted from the experimental data. 

fTabie III contains the resulis caleulated for air "ith a 
¥e 2.0, Sample ecleulations ean be found in Appendix &. 

& comparison of Tables II end III, shoes the correlation 
between the water table anelogy and air flow is within the 


limits of experimental accuracy. 
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SURMARY 


The theory of the enaiogy of two dimensional air flow 





to the fiow of water vith a free surface hes been presented, 

The details of construction of the wetter table, which 
utilizes the water enclegy in the study of air fiow, has 
been described, 

The operation end am application of the tater table 
heve been demonstrated, 

The water table ceuid be employed to a decideG ad- 
vantage by the organizstions which use the few supersonic 
wind tumnels in existance im this coumtry to study the 
flow of high velocity sir. Buch of the work dome at present 
in the wine tumnele could be seeom lished om a water table 
with the resulicnt considerable seving in time ond expense. 

The emphasis in this work has been on supersonic sir 
flow, however, the water table can be used for subsonic 
air fice investigations. It is a fine leboretory apperatus 
fer demonstrating hydraulic and eir floe theory vitueliy. 

The operation of the present seter table noulé be 
enhanced with the use of a depth gauge capable of me:curing 
to 1/1000 ineh, A further refinement in control of mater 
flow is mecéssury. It is suggested that the present sump 
end motor be replaced by «2 coentroliebie direct current 
electric drive virieble speec centrifugel pump, It is sise 
suggested thet inorder to aecursté€iy messure shock sngles 
om photographs, a parellel Light syatem be constructed, 
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SDARSLE CALCULATIONS 


Ren I 

A = Cross sectional srea of water fiow, eyual to he, rt 
H. = Bech number of exter flow, v_/ve 
Vv, = Meter velocity, fps 
Vi = Velocity of propogation of « surface wave, fps 

w= width of weeter flow, ft. 

© + shock engle, dog. 

Jater Flow 

V, = G/A = be = 1.585 ft/sec 
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p/P? = (h./hy)*= (1.987)" = 1.775 
2h, fh, x by fh, = 2085/1.687 2 1.758 


M. = 2% ~ 2) 2\/2¢2.753) = 1.880 
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TabLe JI 
Fifteen Degree fecage 


Experimental Meesurements 
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Fifteen Degree VYedge 


Calcyletec Values 











Table III 
Fifteen Degree Wedge 


Caleulated Values for sir vith Ye ud 
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Notation For Continuity E QuaTion 
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Fig. ll = WATER TABLE (Top View) 
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Abrupt Energy Transformation 


in Flowing Gases 


By N. P. BAILEY GOs ee 


After setting up the basic equations of energy, flow, and 
acceleration, this paper compares theoretical and actual 
plane compression shocks in tubes. By making the single- 
plane-angle shock a special case of a plane shock, curves 
for solving numerical cases are presented. Even though 
most angle shocks encountered in engineering are three 
dimensional, it is felt that the simple theory of single- 
plane shocks is helpful in understanding such phenomena 
without the complications of three-dimensional theory. 
Traverses of three-dimensional angle shocks in nozzle and 
orifice discharges are presented and discussed, and tem- 
perature traverses through shock diamonds in high- 
temperature streams are shown. Although no complete 
explanation exists for thermocouple readings in excess of 
the total gas temperature in such high-temperature 
shocks, their existence is established. As a last case of 
abrupt energy transformation in gas streams, tests of 
combustion with flow in a constant-area tube are given. 
Experiment and analysis agree that the assumptions of 
constant-area flow and steady or continuous flow do not 
appear to be simultaneously tenable when combustion is 
present. 


A 


The following nomenclature is used in the paper: 


NOMENCLATURE 


area, sq ft 
specific heat at constant pressure, Btu per lb per deg F 
specific heat at constant volume, Btu per lb per deg F 
constant of integration 

constants of integration 

total differential 

partial differential 

friction force, lb per ft of length 

acceleration of gravity; 32.2 fps per sec 

rate of heat release, Btu per sec per cu ft 

778.26 ft-lb per Btu 
mechanical work, Btu per lb 
Mach number 

pressure, psfa 

heat added, Btu per lb 

gas constant; 53.3 for air 
temperature, deg Rankine (R) 
time, sec 

specific volume cu ft per lb 
velocity, fps 

weight flow, lb per sec 
distance, ft 

deflection angle, deg 

= ratio of gas specific heats; 


1.395 for cold air 





1 Head, Mechanical Engineering Department, Rensselaer Poly- 
technic Institute. Mem. A.S.M.E. 

Contributed by the Research Committee on Fluid Meters, for 
presentation at the Annual Meeting, New York, N. Y., December 
2-6, 1946, of THE AMERICAN Society OF MECHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. 
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shock angle, deg 
= mass density, slugs per cu ft 


INTRODUCTION 


Mechanical-engineering literature is rich in analytical and 
experimental information on the continuous, steady flow of gases 
but there is need for a wider familiarity with gas flow when 
such discontinuities as compression shocks and combustion with 
flow are present. .The aim of this paper is to present some 
analytical and experimental information on such discontinuous 
and intermittent phenomena. 


FLrow MoMENTUM AND ENERGY CONCEPTS 


In any case of flow, such as illustrated in Fig. 1, the difference 
in the rate at which mass enters and leaves the volume, Adz 
must be balanced instantaneously by the rate of storage, or 


Op 
= —Adz — 
ay 


O(pvA) 
—— d 
Ox : 


At any instant the mass of fluid in the volume (Adz) is (pAdz) 
and, if it receives a change in velocity 


Ov Ov 
dv = —d — dt. . ..ssas eee 2 
% Ox tals ot [2] 
the acceleration is 
dy _ dude , dvdi _ odo , dv “4 
di oxdt § dtdt oe oe ° 


For a wall friction or obstruction force, Fdz, a summation of 
forces in the flow direction gives 
vou 4 ov 
Ox ot 


A useful relationship may be had by multiplying Equation 
[1] by v and adding it to Equation (4) to give 


ov Ov Avodp 
—(pvA) = — pA (2) —Ff— -— [5] 


AoP 
— — dr = Fdx — pAdz 
Ox : 


AP 
Ox 


vo [pA] a 
Ox 7 





PoA 
By adding ee to both sides of Equation [5] and remember- 


OA 
ing that ies 0 for a rigid passage, Equation [5] may be written 
as 


2 
O[PA + pAv?| i O(pAv) i 


PdA —Fdz = a = 


O(pAv) 


For the usual steady-flow case where = 0, Equation 


. ot 
[6] says that the net wall force on the gas in the direction of flow 
is 


PdA — Fdz =d|PA + pAv?]..........05. 


Discussion on this paper will be accepted until February 3, 1947 


2 TRANSACTIONS OF THE A.S.M.E. 


where [PA + pAv?] 1s the total momentum per second in pounds 
passing any section. 


1 
Since |. = — and the Mach number | M = steed 
gRT ygRT 


it is often convenient to express Equation [7]. 


Net wall reaction = PdA — Fdx = d [PA(1 + yM?)].. [8] 





This statement is true only so long as the flow is steady and 
it does not hold when the mass flow (pAv) changes with time. 
However, since any such flow variation must go through a re- 
peated cycle, the average valuc of £ d(pAv) must be zero and 
Equation [7] may be used to evaluate the average wall force even 
when the flow is unsteady or intermittent. This point will be 
covered by experimental work later. 

As gas flows along a channel as in Fig. 1, any thermal energy 





Krew) 


(dQ) that is released or added, and any work done on it (dL), as 
in @ compressor impellér, must be accounted for. Some of it 
will go to change the internal energy (C,dT) and part to do flow 


Vdl 
k _— 
wor ( 7 


a pressure change (dP), a compressible gas will undergo a volume 
PdV 
jf 


and an increase in velocity would induce a kinetic-energy change 


1 
(+. vo) 


This results in an energy-balance equation 


2 
) as the gas flows through a pressure change. With 


change (dV), and do an amount of expansion work 


1 1 1 
dL + dQ =C,dT + — VdP + — PdV + — vdv .... [9] 
J J ag 
A reversible process is characterized by dQ = 0 and 
1 
CydT + Z ati). oc ey es See [10] 


Such a process is thus defined as one where the only internal 
energy change is the inevitable one resulting from compression 
or expansion work. 

A useful form of the energy equation for gencral applications is 


¥ R iI 
————— - dT’ — ydlv...... 1 
G1 7% ae ydv [11] 


dL + dQ = 
For the case of constant total energy (dL = 0 and dQ = 0), if 
Equation [11] is integrated between any velocity v, and corre- 
sponding static temperature 7’, and a final velocity of zero where 
the static temperature is the same as the total temperature 7), 
the result is 





ygRhT + 


Pal 
(y ; ee [12] 


When this is combined with the Mach number definition 


v2 
Sieg sae [13] 
The result is 
Bs i ] 
T (static) = ’ oa eis eave 3) [14] 
1 ee 
2 
For the case of steady flow 
WY = pgAv... 7... [15] 


and this combined with the gas equation and with Equations 
[13] and [14] gives 


WVT) - as (yaa 
ae M R 1+ mu: | os a [16] 


Equation [16] is very useful in all cases of steady flow at con- 
stant total energy. 





PLANE SHOCKS 


For the simple case of steady flow at constant area and con- 
stant total energy with no appreciable wall friction, Equation [8] 
becomes 


PCL + yM?) = Pil + yM,2) = P2(1 + yM2?) = const. . [17] 
and Equation [16] may be written 


Sue —— 
pu sft + 2D we = pa aft + C= 


—l 
1+ oe M,.? = const... [18] 


= P.M, 
When (P2/P;) is eliminated between Equations [17] and [18], 
one obvious solution for the equation of constant area and con- 
stant total energy flow without friction is for (M;, = Mz), but 
when this root (M: — M)) is factored out, the result is 


_ Gaaaw M,? 
M,? = —— es [19] 
—] 
yM,? 2, ) 
and 
27M (7— 
P,/P, = ——————, 2 
2/ Py G41) [20] 


This means that for any superacoustic Mach number M,, 
there is a subacoustic Mach number Mg, which satifies the condi- 
tions of flow. This defines a plane compression shock which 
is a discontinuity, occurring theoretically in an extremely short 
distance. Whether or not it will occur in any case depends 
upon operation conditions. Since channel friction? causes 
superacoustic flow to approach the acoustic, for any initial 
Mach number there is a maximum flow distance for which shock- 
less flow is possible. 

A shock may also be induced in a shorter channel by applying 
the correct back pressure. For any initial Mach number and 
channel friction, each position of the shock uniquely determines 
one back pressure. This is illustrated in Fig. 2, where A —1—B 
represents flow that is all superacoustic. If the back pressure 





2 “The Thermodynamics of Air at High Velocities,” by N. P. 
Bailey, Journal of the Aeronautical Sciences, vol. 11, July, 1944, p. 
Lat. 
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s raised somewhat, the flow adjusts to it by a plane compres- 
sion shock at some point 1, going to subacoustic flow at 2, and 
discharging subacoustically at C. Equation [19] indicates that 
Oey) 

27 
very large, and it cannot go below that value. 

The magnitudes of the pressure rise and Mach number change 
of actual compression shocks check well? with theoretical values 
from Equations [19] and [20]. However, instead of being com- 
pleted in a negligible distance, the pressure rise extends over a 
distance of from one to five pipe diameters. 

This is illustrated by the test data shown in Fig. 3. Both wall 
and center pressure traverses were taken for the conditions shown 


the Mach number after shock approaches as M, becomes 
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Fic. 3 CoMPprESSION SHOCK IN A TUBE 


Pov 


and after a series of angle shocks in the first 1.15 in. of the tube, 
the wall and center pressures became the same, and a plane com- 
pression shock was initiated. ‘The Mach number as calculated 
from Equation [16] for an initial pressure of 14.33 in. of mercury 
was M, = 1.37. Similarly at the final pressure of 29.50 in., 
M, was 0.75. From Equation [19], a shock, started at M; = 
1.37, should have ended at Mz = 0.748, and from Equation [20] 
the pressure should be P, = 29.50. 

The only discrepancy lies in the fact that the pressure rise is 
distributed over a length of 1.3 in. instead of being an abrupt 


change as illustrated in Fig. 2. The explanation of this lies in 


the effect of boundary layer and the use of a manometer to 
measure pressure. Since the flow in the boundary layer is sub- 
acoustic, no compression shock and corresponding pressure rise 
can occur in it. This means that at any instant the large pres- 
sure rise across the shock front is short-circuited through the 
boundary layer. This immediately causes the shock to col- 
lapse and form at some other point, only to collapse again. 
Since the manometer reading is a time average of the pressure 
at any point, a value of one half the pressure rise, occurring at 
1.4 in., Fig. 3, merely means that the compression shock is 
upstream from that point one half the time and downstream 
from there the other half of the time. The shock is constantly 
dancing back and forth in the tube, never being further upstream 
than the 1.17-in. point and never further downstream than 2.5 
in. 

This explanation is borne out by high-speed photography, as 
well as by the fact that an open-end impact tube tunes to a shrill 
whistle in such a shock region. Consequently, the simple pic- 
ture of Fig. 2 will have to be modified to show the compression 
shock (1-2) as sweeping back and forth in the tube through a 
distance that varies with boundary-layer thickness, tube size, 
and initial Mach number. 


ANGLE SHOCKS 


When supersonic gas flow encounters an obstruction which 
calls for it to change its flow direction, an angle shock will occur. 
This is because the approach velocity is so high that the first 
molecules to find themselves in trouble are not capable of sending 
a distress signal upstream. This means that each successive layer 
of gas molecules must run into the same trouble with no warning. 

This is illustrated in Fig. 4 for the simple two-dimensional 
case where supersonic parallel flow at a velocity m, must turn 
through an angle a, and flow parallel in a new direction. This 
entire change takes place in a shock front at an angle @, from the 
original flow direction. Such an angle shock may be made into a 
special case of the plane shock, discussed previously, by viewing 
it as a pronmias velocity (v; sin 6), going through a plane shock to a 
velocity v2’, and this all occurring in a flow field having a uniform 
velocity (v, cos 9) parallel to the shock front. This can be. ‘done 





Trig. 4 


4 TRANSACTIONS OF THE A.S.M.E. 


beeause there is a pressure rise normal to the shoek front but 
no pressure gradient to give a veloeity ehange parallel to the 
shoek front. 

From Fig. 4 the condition that the normal component of the 
velocity must be reduced to a value v9’, that will eause the final 
velocity vz, to be turned through an angle a, may be stated as 


Dv» ; 


- 


——- = tan (0— a) See. ae [21] 
v, eos 6 


From Equations [13] and [14] 


V VGRP o1 M, cos 6 








v, cos 8 = Se Mees dutta [22] 
all 
\' i vet) M,? 
2 
Also, sinee v;’ = », sin 6, for the same static temperature 7, = 7, 
Equation [13] says that 
M,’ = Misin @........ ae [23] 
From Equation [19] then 
-_ —1 
1+ a M,’2 1+ u M,?sin20 
OO .. [24] 
, ope : (y a) 
y M,’?— zs a 


Before (v2’) can be evaluated from Equations [22] and [24], the 
total temperature 7’. = 1,’ must be expressed. The static 
temperature 7,’ = T; may be evaluated from Equation [14] as 


To 

T, = T,' = Ba (25 ] 
(ya) 
1+ M,? 
2 
Since, from Equation [22] 
RT... M 6 

v,' — Uy Sie MV vgs Mi sine = RE senate [26] 
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and from Equation [12] the total temperature Jo,’ is 


(ial a 
a ‘= ne = Tic _———— 2a 
l CS QvgR v1 [27] 






This gives the total temperature 752’ as 


Ta (y a 1) Ls M, 2 sin? 6 1 ; 
_ = te aC 
1+ 2= 2 wy: 2a hn 


Using Equations [24] and [28° in 


/ 
fm a 


RT o' Mo’ | 
Vo! & _V gh Pon! Met 11. en (29 


gives 





(y= ; 
—— 1 —— M,? 26 
, V ygRT 1 | + 2 Lei 


oo  S= SSSSSSS—eeee Oe 


\' i ee 2 1) M,? (y+) M, sin 0 


From Equations [30] and [22] in [21] 


—] 
[ + wSs 5 ) M,? sin? | 


fin (@ =="0) & 
an (0 — a) 80" ee 


Equation [31] defines the shock angle @ for gas flow at an initia 
Mach number M,, being deflected through an angle a. 

Using M,’ = M, sin @ in Kquation [20] gives the statie-pres 
sure-ratio rise through the single-plane-angle shock as 


(y= 
(ya 


The resultant final veloeity v; is given by 


i, 2 7 M,?sin?6 
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Using Equation [21] for v2’ and Equation [22] for », cos @ in 
Equation [33] 





re \/ hae M, cos 6 | [34] 
— 1) 
-——§ M,? cos (6 — a) 
Z 
From Equations [12] and [13] 
1 
M, = Sent REC ECn > >: 135] 
y9RT 2  (y—)) 
Vo? 2 


Since 2 and v2 are total velocities and a compression shock is at 
constant total energy, Z'n = Z'e 


i 
M: TTT Tn nee een en Oe |. is [36] 
=, , 
1+ SD ay:e cos? (@ — a) 
2 mo!) 
M2 cos? @ 2 


This gives the final Mach number Mb, after an angle shock from 
an initial value of My, when a le-plarie flow is deflected through 
an angle a. 

Fig. 5 shows values of shock angle 6, plotted against the deflec- 
tion angle a, for various values of initial Mach number as calcu- 
lated from Equation [81]. For each value of M, there is a 
maximum deflection angle for which angle shock conditions can 
be satisfied and for any increased deflection angle a plane shock 
will result. This maximum point represents the conditions 
where a wedge in a free stream would cease to have angle shocks 
from its nose and would have a bow shock out in front of it. 

Fig. 6 shows corresponding values of shock-pressure ratio for 
var. 4s conditions, anu Fig. 7 shows the final Mach number after 
an angle shock. Whereas the flow after a plane shock is always 
subacoustic, it may be either above or below sonic in the case of 
angte shocks. 

When supersonic air flows past a wedge, the angle shocks are 
stably anchored at the nose of the wedge, but when a boundary 
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layer is present, th» shock pressure rise is distributed over an 
appreciable distance as in the case’of a plane shock. 

In Fig. 3, where the first angle shock originated at the wall, the 
pressure rise extended from 0.14 in. to 0.22 in. The first pres- 
sure rise indicated by the center trav- 
erse spread between 0.56 in. and 0.72 
in. In the next angle shock at the wall 
where the boundary layer was thicker, 
the pressure rise extended from 0.80 
in. to 0.95 in. This would indicate that 
angle shocks emanating from a channel 
wall are unstably anchored to a sub- 
sonic boundary layer the same as plane 
shocks. 


SHocks at High TEMPERATURES 


When a gas flows through an orifice or 
a parallel-walled nozzle under a pressure 
ratio greater than that needed to pro- 
duce sonic velocity, the sudden adjust- 
ment of the gas pressure to the lower 
back pressure produces a variety of angle 
shock patterns. Fig. 8 shows static and 
impact pressure center traverses for 
such flow from a simple nozzle, to- 
gether with the probable shock and 
rare-faction wave pattern. 

The fact that each pressure rise is 
spread over a distance of 0.15 in. indi- 
cates that, owing to the boundary layer 
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in the nozzle, the shock pattern was not stable and stationary 
in space. Further evidence of such instability was furnished by 
the whistling of the impact tube in the region from 0.05 in. 
inside the nozzle to 0.10 in. outside of it. 

The floss of impact pressure and the subsequent recovery of it 
between 0.15 and 0.40 in. marks the region where the flow was 
sufficiently supersonic to produce a serious shock in front of 
the impaet tube. The distance of this shock bow in front of the 
impact tube appears to have been about 0.08 in. 

It is usually assumed that the first angle compression shock 
emanates from the nozzle wall, but the traverse of the jet from 
an orifice, shown in Fig. 9, indicates that angle compression 
shocks can originate from a jet boundary. Sonie velocity was 
not reached until 0.2 in. from the orifice and no disturbance 
reached the center traversing tube in the first 0.33 in. At that 
point an angle shock and the subsequent rarefaction were fol- 
lowed by a plane or bow shock in front of the tube shown in the 
stream. 

Such shock patterns as shown by the curves in Figs. 8 and 9, 
respectively, are quite common, and alternate compressions and 
rarefactions usually persist for 6 or more oscillations before they 
are damped out. With low-tempcrature gases, they are detecta- 
ble only by making pressure traverses or by schlerin photo- 
graphs. 

When these shock patterns occur in high-temperature prod- 
ucts of combustion, the high-pressure-shock diamonds are directly 
visible and are usually a beautiful blue color. Such a view is 
illustrated in Fig. 10. Fig. 11 shows test results for thermo- 
couple readings taken in such a shock region in a stream of 
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products from an orifice approximately 2 in. diam. The first 
reaction to such results as shown in Fig. 11 is that of utter dis- 
belief. However, the chance of unusual instrumentation errors 
has been systematically eliminated by repeating such test runs 
under a variety of conditions, always with the same results. 

The next idea is that there is unburned fuel being discharged 
from the orifice, which burns on the surface of the thermocoup 
and produces a local temperature rise. The chance of this being) 
the case is very small] for two reasons: (a) The original burnél 
temperature checks a heat balance on the air and fuel supplhie 
within 2 per cent; (b) the appearance and color of the shoe 
diamonds are not at all altered by the insertion of a thermocouple 

Another possible explanation is that some unstable oxides « 
nitrogen are formed in the region of shock where there are rate 
of deceleration greater than 2,000,000 g. This again has bee 
eliminated by burning propane with pure oxygen and observing! 
exactly the same phenomenon. 

Fig. 12 indicates that it is a phenomenon definitely associate 
with shocks. The data for the curves in Figs. 12 and 13, re 
spectively, were taken in a 3/,¢.-in. jet with a fairly large thermo= 
couple (1/s-in. sheath), so the cooling effect was quite pronouncec 
As a result, the excess temperatures shown are much smalle 
than those shown for the large jet in Fig. 10, but they are mor 
analytical. 

Fig. 12 indicates further that this phenomenon is definitel, 
geared to the strength of the compression shock, for it com= 
pletely disappears at orifice or burner pressure ratios where’ 
shocks are not present. It would naturally be expected that an_ 
unshielded thermocouple in a small stream of hot gases woul 
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‘iq. 10 HiGH-PRESSURE-SHOCK DIAMONDS IN H1GH-TEMPERATURE 
Propuctrs oF COMBUSTION 
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Fic. 12 (Above) THERMOcoUPLE READINGS IN SHock DiaMonps 
Usinc CompustTion Propucts or PROPANE AND AIR: EFFECT OF 
BURNER PRESSURE RaTIO ON Excess TEMPERATURE READINGS 


THERMOCOUPLE READINGS IN SHOCK DIAMONDS 
ISFFECT OF 
READINGS, 


Fic. 13 (Below) 

Usine COMBUSTION PRODUCERS OF PROPANE AND AIR: 

BuRNER TEMPERATURE ON Excess TEMPERATURE 
BURNER PRESSURE Ratio = 3.62 


read from 100 to 200 deg F low owing to radiation losses and 
velocity errors, and that is what happened below pressure ratios 
of 1.9. 

Fig. 13 indicates that vain the limits of accuracy of the test 
data, the amount of excess temperature is fairly independent of 
burner temperature above 1100 F. No test data are available 
for initial temperatures between 1100 F and 80 F. However, 
repeated attempts to find excess thermocouple readings in 
shocks at a temperature at 80 F failed. 

Since the thermocouple metal reaches a temperature greater 
than the total temperature of 
the gas before it is in energy 
equilibrium, it is evident that 
it is being bombarded by some- 
thing besides gas molecules. 
The violent deceleration of the 
hot gases in the shock could 
produce such a high concen- 
tration of ions that their bom- 
bardment and _ neutralization 
at the thermocouple metal 
surface could cause the metal 
to reach a very high tempera- 
ture before energy equilibrium 
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is reached. This phenomenon deserves a thorough investiga- 


tion. 
CoMBUSTION AND FLow 


Combustion of air-fuel mixtures is customarily done with no 
appreciable flow velocity and too little is known of the mechanism 
involved when combustion oeeurs in a high-velocity air stream. 
The most commonly discussed case of combined combustion and 
flow is that of burning in a constant-area, passage which is such a 
short distance that wall friction forces may be neglected. 

If it is also assumed that the area is constant and the flow is 
steady, Equation [7] gives the condition for constant momentum 
per second as 


Ls + pve” airciac dy circiae fs ureiroth nemo oas [37] 


solving for the increase in veloeity pressure 


Pi, + piv? = 


This says that the gain in velocity pressure is only one half the 
drop in statie pressure. It was to find out, if possible, how the 
gas stream goes about losing the rest of this statie pressure with- 
out wall friction that the following tests were made: 

The 1%/s-in-diam thin-walled tube, used in Figs. 14 and 15, 
was lined up axially in the air stream from a 3-in-diam nozzle. 
Hydrogen for heating was introduced through the 3/s-in-diam 
tube shown. By regulating the size and arrangement of hy- 





3 Such a program is being carried on by the Mechanical Engi- 
neering and Physics Departments of the Rensselaer Poly technic 
Institute. 
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bic.14. CoMBUSTIONINSTRAIGUT TUBE; I¢xpLOSIVE BURNING 


drogen jets and the hydrogen pressure, a uniform temperatun 
could be reached aeross the entire cross section. Statie anc 
total pressures was obtained with !/jo-in. traversing tubes, ant 
temperatures with an unshielded thermocouple of the sam 
maximum diameter. All values shown represent averages 0 
the annulus around the hydrogen tube. 

In general, two quite different sets of conditions could be ob 
tained. When the heating was uniform across the tube the re 
sults in Fig. 14 are representative of several tests. The burnir 
was a series of explosions at the resonant frequency of the tub 
This frequeney was high enough for manometers and therm 
couples to give steady readings, but it was a violently noisy forn 


of combustion. The thing which characterizes these result: 
vp 


shown in Fig. 14 is that the velocity pressure (= inereases 1 


the region where the temperature is rising. 

However, when the hydrogen tube was placed slightly off cente 
and the tube wall became a bit hotter on one side, the eombustio: 
beeame anchored and was very quiet and orderly. The traverse 
in Fig. 15 are typical of these tests, and they are characterize 
by the fact that practieally all of the velocity-pressure inereasé 
occurred ahead of the temperature rise. The chief loss in tota 
pressure occurred in the burning region where the velocity pres 
sure remained essentially constant. 

This would indicate some form of jet separation in the tub 
ahead of the combustion, followed by combustion at essentially 
constant velocity pressure as the flow again filled the tube 
Fig. 16 shows the velocity plotted against absolute temperature 
and for the tube to flow full, the velocity would be essentially 
proportional to temperature, that is, if the weight of hydrogen 
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Fic. 15 Compustion IN StrAicuT Tune; Quiet BURNING 


idded and the small pressure change are neglected. The points 
yefore and after combustion do lie very close to such a line, but 
she velocity in the burning region is much higher. This furnishes 
. good indication that the tube was not flowing full. 

In the second section of the paper it was concluded that even 
vith unsteady flow, Equation [383] would be true for averaged 
‘eadings. For the case of the curve in Fig. 14, for intermittent 
Jow and explosive burning, the ratio of velocity-pressure increase 
50 static-pressure decrease was 0.65 instead of the theoretical 0.50. 
{n the steady-flow case in Fig. 15 it was 0.57, and for another 
series of eight tests it varied from 0.52 to 0.58. Since Equation 
{38] does not account for the addition of hydrogen to the flowing 
zas stream this is not a bad check. 

It is easy to assume that just because combustion occurs in a 
constant-area tube the burning occurs at constant area. The 
imited amount of test work at hand would seem to cast doubt 
yn the validity of the simultaneous assumptions of steady flow 
and constant-area combustion. 

Some of the implications of these assumptions may be had by 
assuming a rate of heat release of H Btu per sec per cu ft of vol- 
ume, and stating that for steady operation an amount of energy 
AHdz must be transported out of the volume Adz, each second; 
or, from the energy equation 







1 


W 
HAdz = WC,dT + ce BHU nee onan [39] 
For acceleration, with no wall friction 
1 W 
AdP = —pAdz sage —pAvdy = ——dv...... [40] 
dt g 
For steady flow 
Pe 
= eet ee ret 4] 
W = pAgu RT Av [41] 


Integrating Equation [40] 
Ww 
AP =——-v+(Q Seta a eee Reka ss) ele l stats [42] 
g 
For any initial condition P; and 1, C, may be evaluated to give 
A 
a a Ee ile. [43] 


From Equations [40] and [43] 





gA Ag 
=— — == Pees 44 
vdy wor E + a (P; Py | [44] 
Differentiating Equation [41] 
A A A 
oo —— i pee Een oso os 
RW d(Pv) pwr + Rw Pdv [45] 
From Equations [40], [43], and [45] 
A Ag gA*PdP 
citM = rw ee E + WV (P,— P) i RW [ 46] 
From Equations [44] and [46] in [89]. 
ie HW 
meno 8~OCi\lwn.~SC—sSC [47] 
oy = a: 
a i(¢ aq ail | R Pt 
Since, from Equation [41] 
PiA 
Wu = Ww, = ai Te [48] 
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Fie. 16 VELocITY-TEMPERATURE TRAVERSE FOR CuRVE, Fic. 15 


and from definition 





= Bo aw gan Se 49 
ygRT, a 
and 
Cine. 1 
ee eee 50 
R J EB J(y¥—1) (50) 


Equations [48], [49], and [50] with Equation [47] given 


dP HJIW(y — 1) 1h 
a ee se .s a [ow 
Jn.) PygA ae 


7? 
(y + 1) pr Me + 1) 
1 


This type of a pressure gradient has all the earmarks of in- 
stability, for (dP/dx) becomes a greater negative value as P be- 
comes smaller. This would indicate that the effect would tend 
to cumulate until (dP/dx) becomes minus infinity which would 
be a discontinuous front. Consequently, analysis appears to 
bear out the experimental conclusion that true constant-area 
burning cannot also be a steady-flow phenomenon. 
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